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Full LEP-I data collected by the ALEPH detector during 1991-1995 running are analyzed in order to measure 
the r decay branching fractions. The analysis follows the global method used in the published study based on 
1991-1993 data, with several improvements, especially concerning the treatment of photons and 7r°’s. Extensive 
systematic studies are performed, in order to match the large statistics of the data sample corresponding to 
327148 measured and identified r decays. Preliminary values for the branching fractions are obtained for the 2 
leptonic channels and 11 hadronic channels defined by their respective numbers of charged particles and 7r°’s. 
Using previously published ALEPH results on final states with charged and neutral kaons, corrections are applied 
so that branching ratios for exclusive final states without kaons are derived. Some physics implications of the 
results are given, in particular concerning universality in the leptonic charged weak current, isospin invariance in 
ai decays, and the separation of vector and axial-vector components of the total hadronic rate. 


1. Introduction 

A complete and final analysis of r decays is pre¬ 
sented using a global method. All data recorded 
at LEP-I with the ALEPH detector are used, 
thus providing an update of those previous results 
which were based on partial data sets. The in¬ 
crease in statistics —the full sample corresponds 
to about 2.5 times the luminosity used in the last 
published global analysis |l|,|| not only allows 
for a reduction of the dominant statistical error 
but, more importantly, provides a way to bet¬ 
ter study possible systematic biases and to even¬ 
tually correct for them. Several improvements 
of the method have been introduced in order to 
achieve a better control over the most relevant 
systematic uncertainties: simulation-independent 
measurement of the rr selection efficiency, im¬ 
proved photon identification especially at low en¬ 
ergy where the separation between photons from 
r decays and fake photons from fluctuations in 
hadronic or electromagnetic showers is delicate, 
a new method to correct the Monte Carlo sim¬ 
ulation for the rate of fake photons, and stricter 
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criteria for channels with low branching fractions. 
For consistency and in order to maximally profit 
from the improved analysis all data sets recorded 
from 1991 to 1995 have been reprocessed. The re¬ 
sults presented in this paper thus supersede those 
already published in Ref. [|||TJ|^]. Only the mea¬ 
surements on final states containing kaons, which 
were already based on the full statistics, remain 
unchanged 


2. Experimental method 

2.1. The data and simulated samples 

A detailed description of the ALEPH detector 
can be found elsewhere 10 Il| . 

Tau-pair events are simulated by means of a 
Monte Carlo program which includes initial state 
radiation computed up to order a 2 and exponen¬ 
tiated, and final state radiative corrections to or¬ 
der a ||. The simulation of the subsequent r 
decays also includes single photon radiation for 
the decays with up to three hadrons in the final 
state. The longitudinal spin correlation is taken 
into account |[3]]. This simulation, with the de¬ 
tector acceptance and resolution effects, is used 
to initially evaluate the corresponding relative ef¬ 
ficiencies and backgrounds. It also includes the 
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tracking, the secondary interactions of hadrons, 
bremsstrahlung and conversions. For all these ef¬ 
fects, detailed comparisons with relevant data dis¬ 
tributions are performed and corrections to the 
MC-determined efficiencies are derived. 

The data used in this analysis have been 
recorded at LEP I in 1991-1995. The numbers 
of detected r decays are correspondingly 132316 
in 1991-1993 and 194832 in 1994-1995, for a to¬ 
tal of about 3.3 • 10 5 . The ratios between Monte 
Carlo and data statistics are 7.3 and 9.7 for 1991- 
1993 and 1994-1995 periods, respectively. Monte 
Carlo samples are generated for each year of data- 
taking in order to follow as closely as possible the 
status of the detector components. 

2.2. Selection of rr events 

The principal characteristics of rr events in 
e + e _ annihilation are low multiplicity, back-to- 
back topology and missing energy. Each event 
is divided into two hemispheres by an energy 
flow algorithm JTTJ which calculates all the vis¬ 
ible energy avoiding double-counting between the 
TPC and the calorimeter information. The jet 
in a given hemisphere is defined by summing 
all the four-momenta of all energy flow objects 
(charged and neutral). The energies in the two 
hemispheres including the energies of photons 
from final state radiation, E\ and E 2 , are use¬ 
ful variables for separating Bhabha, /i/x and 77 - 
induced events from the rr sample, while the rel¬ 
atively larger jet masses, wider opening angles, 
and higher multiplicities indicate Z —> qq events. 

All these features are incorporated in a stan¬ 
dard selector used extensively in ALEPH |p^,p|2]. 
In the previous analyses [£j]j|] additional cuts had 
been introduced in order to further reduce the 
contamination from ee and /r/i processes. In the 
present work it was chosen to simplify the pro¬ 
cedure in order to conveniently measure selec¬ 
tion efficiencies on the data, at the expense of a 
slightly larger background contamination which 
is anyway also measured in the data sample as 
explained later. 

We use the ’break-mix’ method introduced for 
the determination of the rr cross section [Q to 
measure the efficiency of all the selection cuts. 
For every cut, one hemisphere of the event is 


chosen judiciously so that it is unbiased with re¬ 
spect to the cut under study and free of non-r 
backgrounds. This procedure selects the opposite 
hemisphere as an unbiased r decay which is then 
stored away. Pairs of selected hemispheres are 
combined to construct a rr event sample built 
completely from data. This sample is used to 
measure the efficiency of the given cut. 

The measured efficiencies are found to be very 
close to those obtained by the simulation, devi¬ 
ations being at most at the few per mille level. 
This situation stems from the facts that the r 
decay dynamics is —apart from small branching 
ratio channels— very well known, the selection 
efficiencies are large and the simulation of the de¬ 
tector is adequate. The overall selection efficiency 
of rr events is 78.8 %. This value increases to 
91.9 % when the rr angular distribution is re¬ 
stricted to the detector polar acceptance, giving 
a better indication for the efficiency of the cuts 
designed to exclude non-rr backgrounds. In ad¬ 
dition, when expressed relatively to each r decay, 
the selection efficiencies are weakly dependent on 
the final state, with a total relative span of only 
10 % for the 13 considered decay topologies. 


2.3. Estimation of non-rr backgrounds 

A new method —already implemented for the 
measurement of the r polarization Jl5[ - has 
been developed to directly measure in the final 
data samples the contributions from the major 
non-r backgrounds: Bhabhas, /i + /i~ pairs, and 
77 —> e + e - ,/r + /i“, and hadrons events. The pro¬ 
cedure does not require an absolute normalization 
from the Monte Carlo simulation of these chan¬ 
nels, only a qualitative description of the distri¬ 
bution of the discriminating variables. The basic 
idea is to apply cuts on the data in order to reduce 
as much as possible the rr population while keep¬ 
ing a high efficiency for the background source 
under study, i.e. the reverse of what is done in 
the rr selection. 

The non-r backgrounds in each channel are 
listed in Table [|. They amount to a total fraction 
of (1.23 ± 0.04) % in the full data sample. 
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2.4. Charged particle identification 

A ’good’ track is defined to have a momentum 
greater than or equal to 0.10 GeV/c , | cos(0)| < 
0.95, at least 4 hits in the TPC, and its min¬ 
imum distance to the interaction point within 
2 cm transversally and 10 cm along the beams. 
In classifying r decays, only good tracks are used, 
after removing those identified as electrons which 
are used to reconstruct converted photons; the 
electrons identified as bad tracks are also included 
in the reconstruction of conversions. 

Charged particle identification is achieved with 
a likelihood method incorporating the informa¬ 
tion from the relevant detectors. In this way, each 
charged particle is assigned a set of probabilities 
from which a particle type is chosen. No attempt 
is made in this analysis to separate kaons from 
pions in the hadron sample since final states con¬ 
taining kaons have been previously studied . 

Eight discriminating variables are used in the 
identification procedure: dE/dx in the TPC, 
two estimators (transverse and longitudinal) of 
the shower profile in ECAL, the average shower 
width measured with the HCAL tubes in the fired 
planes, the number of fired planes among the last 
ten, the energy measured with HCAL pads, the 
number of hits in the muon chambers (±4tr-wide 
around the track extrapolation, where a is the 
standard deviation expected from multiple scat¬ 
tering), and finally, the average distance (in units 
of the multiple-scattering standard deviation) of 
the hits from their expected position in the muon 
chambers. 

The performance of the particle identification 
has been studied in detail using control samples of 
Bhabha events, ///i pairs, 77 -induced lepton pairs 
and hadrons from 7 r -tagged r decays over the 
full angular and momentum range |I|] . Measured 
efficiencies and misidentification probabilities are 
given in Figs. and |^. 

2.5. Photon identification 

The high collimation of r decays at LEP ener¬ 
gies quite often makes photon reconstruction dif¬ 
ficult, since these photons are close to one another 
or close to the showers generated by charged 
hadrons. Of particular relevance is the rejec¬ 
tion of fake photons which may occur because 
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Figure 1. Hadron identification efficiency and 
misidentification probabilities obtained from the 
tt Monte Carlo, corrected from data using the 
control samples, for the 1994-1995 data set. 


of hadronic interactions, fluctuations of electro¬ 
magnetic showers, or the overlapping of several 
showers. These problems reach a tolerable level 
thanks to the fine granularity of ECAL, in both 
transverse and longitudinal directions, but nev¬ 
ertheless they require the development of proper 
and reliable methods in order to correctly identify 
photon candidates. 

A cluster in ECAL is accepted as a photon can¬ 
didate if its energy exceeds 300 MeV and if its 
barycentre is at least 2 cm away from the closest 
track extrapolation. 

A likelihood method is used for discriminat¬ 
ing between genuine and fake photons. For every 
cluster, an estimator P 1 is defined, P 7 = 0,1 cor¬ 
responding to fake and real photons,respectively. 
It is constructed using probability densities ob¬ 
tained by simulation, but corrected through de- 
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Figure 2. Lepton identification efficiencies ob¬ 
tained from the tt Monte Carlo, corrected from 
data using the control samples, for the 1994-1995 
data set. 






tailed comparisons between data and Monte 
Carlo. 

Discriminating variables for each photon can¬ 
didate used are the distance to the charged track, 
the distance to the nearest photon, a parameter 
from the clustering process, the fractions of en¬ 
ergy deposition in ECAL stacks 1, 2 and 3, and a 
parameter related to the transverse energy distri¬ 
bution. Major improvements were introduced at 
this stage in the analysis compared to the previ¬ 
ous one Q , mainly in the choice of variables and 
in the use of energy-dependent reference distribu¬ 
tions. 

Fig. H shows the comparison of the photon iden¬ 
tification probability (P 7 ) before and after intro¬ 
ducing energy-dependent reference distributions 
and optimized variables. A spectacular improve¬ 
ment is observed in genuine and fake photon dis¬ 
crimination: at low energy, a clear contribution of 
genuine photons can be seen, while at high energy, 
a small, but well identified, fake photon compo¬ 
nent shows up. 

Better photon energy calibration is also 


Figure 3. Comparison of the probability distri¬ 
butions of a cluster to be identified as a genuine 
photon before (histogram) and after (error bars) 
introducing energy-dependent reference distribu¬ 
tions and an improved set of discriminating vari¬ 
ables. 


achieved compared to the previous analysis. The 
procedure, aiming at a relative calibration be¬ 
tween data and simulation, is implemented in sev¬ 
eral steps. First, a calibration is done using elec¬ 
trons from r decays, treating the ECAL infor¬ 
mation as for a neutral cluster. This method is 
reliable only above 2 GeV because of the electron 
track curvature in the magnetic field. Then the 
final calibration is achieved by comparing the re¬ 
constructed 7r° mass and resolution as a function 
of energy. 

Photon conversions are reconstructed following 
the previous analysis ||]. 
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2.6. 7T° reconstruction 

Three different kinds of 7r°’s are reconstructed: 
resolved 7r° from two-photon pairing, unresolved 
7T° from merged clusters, and residual 7 r° from 
the remaining single photons after removing ra¬ 
diative, bremsstrahlung and fake photons with a 
likelihood method. 

A 7T° estimator is defined ||, taking into ac¬ 
count the photon identification probabilities and 
the compatibility of the pair with the 7r° mass 
hypothesis. 

As the 7 r° energy increases it becomes more dif¬ 
ficult to resolve the two photons and the clus¬ 
tering algorithm may yield a single cluster. The 
two-dimensional energy distribution in the plane 
perpendicular to the shower direction is examined 
and energy-weighted moments are computed. As¬ 
suming only two photons are present, the sec¬ 
ond moment provides a measure of the 77 invari¬ 
ant mass pT|. Clusters with mass greater than 
0.1 GeV/c <ire taken as unresolved 7r°’s. 

After the pairing and the cluster moment anal¬ 
ysis, all the remaining photons inside a cone of 
30° around the jet axis are considered. The ra¬ 
diative and bremsstrahlung photons are selected 
using the same method as described in the pre¬ 
vious analysis (Q. Radiative and bremsstrahlung 
photons are not used in the r decay classification 
discussed below. The behaviour of these estima¬ 
tors was studied in Ref. |J. The agreement be¬ 
tween the number of bremsstrahlung photons in 
data and simulation is good, but concerns mainly 
the electron r decay channel where this contri¬ 
bution is important (however it does not affect 
the rate). Bremsstrahlung photons (i.e. radiation 
along the final state charged pions) in hadronic 
channels are at a much lower level and they are 
difficult to pick up unambiguously in the data. To 
estimate the effect of this contribution one there¬ 
fore largely relies on the description of radiation 
at the generator level in the Monte Carlo ||l6|. 

2.7. Decay classification 

Each r decay is classified topologically accord¬ 
ing to the number of hadronic charged tracks, 
their particle identification and the number of 
7 r°’s reconstructed. While for 1-prong and 5- 
prong channels the exact multiplicity is required, 


the track number in 3-prong channels is allowed 
to be 2, 3 or 4, in order to reduce systematic ef¬ 
fects due to tracking and secondary interactions. 

The definition of the leptonic channels requires 
an identified electron or muon with any number of 
photons. Some cuts on the total final state invari¬ 
ant mass are introduced to reduce feedthrough 
from hadronic modes Q .Also some decays with at 
least two good electron tracks are now included, 
when one or more of such tracks are classified as 
converted photons. 

In the previous analysis 0^ the 3 /i27t°, 3/i37t° 
and h4:TT° channels suffered from large back¬ 
grounds and consequently had a low signal-to- 
noise ratio. Most of these backgrounds are due 
to secondary interactions of the hadronic track 
with material in the inner detector part. In order 
to improve the definition of these channels the 
following steps are taken: (1) require the exact 
charged multiplicity n c h = 3 (instead of 2, 3, or 4) 
for 3h2n 0 and 3h3n°, (2) demand a maximum im¬ 
pact parameter of charged tracks less than 0.2 cm 
(instead of 2 cm) for 3h3n°, and (3) require the 
number of resolved 7r°s in 3/i37r° to be 3 or 2, and 
4 or 3 in /i47t°. With these tightened cuts the 
signal-to-noise ratio improves significantly with a 
small loss of efficiency. 

It should be emphasized that all hemispheres 
from the selected rr event sample are classified, 
except for single tracks going into an ECAL crack 
(but for identified muons) or with a momentum 
less than 2 GeV (but for identified electrons and 
hemispheres with at least one reconstructed 7r°). 
These latter decays, in addition to the rejected 
ones in the 3/i27r° or 3h3n° channels are put in a 
special class, labelled 14, which therefore collects 
all the non-selected hemispheres. By definition, 
the sample in class-14 does not correspond to one 
physical r decay mode. In fact, it follows from the 
simulation that this class comprises about 21% 
electron, 27% muon, 41% 1-prong hadronic and 
11% 3-prong t decays. However, consideration of 
class-14 events can test if the rejected fraction is 
correctly understood, as discussed later. 

The numbers of r’s classified in each of the con¬ 
sidered decay channels are listed in Table [[]. 

The KORAL07 generator jlj] in the Monte 
Carlo simulation incorporates all the decay modes 
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Table 1 

Numbers of reconstructed and estimated non-r 
background events in 1991-1993 and 1994-1995 
data sets in the different considered topologies. 


class 

91-93 

94-95 

e 

22405 

598 ± 46 

33100 

745 ± 58 


22235 

409 ± 45 

32145 

380 ± 40 

h 

15126 

93 ±11 

22429 

100 ± 13 

hn° 

32282 

141 ± 22 

49008 

178 ± 26 

h2n° 

12907 

44 ±9 

18317 

81 ± 16 

hAn 0 

2681 

26 ±7 

3411 

35 ±9 

hAn 0 

458 

12 ±3 

499 

19 ±5 

Ah 

11610 

87 ±20 

17315 

129 ± 30 

Ahn° 

6467 

97 ±23 

9734 

165 ± 39 

Ah2n° 

1091 

27 ±7 

1460 

36 ± 10 

AhAn 0 

124 

13 ±4 

150 

25 ±7 

5 h 

60 

3 ± 1 

105 

7 ± 2 

5hn° 

36 

16 ±5 

59 

21 ±6 

14 

4834 

249 ± 38 

7100 

303 ± 52 

sum 

132316 

1815 ±86 

194832 

2224 ± 107 


considered in this analysis, except for the hAn° 
decay channel. In the latter case, a separate gen¬ 
eration was done where one of the produced r 
is made to decay into that mode using a phase 
space model for the hadronic final state, while 
the other r was treated with the standard de¬ 
cay library. The complete behaviour between the 
generated decays and their reconstructed coun¬ 
terparts using the decay classification is embod¬ 
ied in the efficiency matrix. This matrix e 3 , L gives 
the probability of a r decay generated in class 
j to be reconstructed in class i. Obtained ini¬ 
tially using the simulated samples, it is corrected 
for effects where data and simulation can possibly 
differ, such as particle identification, as discussed 
previously, and photon identification as affected 
by the presence of fake photons. 


and bias the efficiency matrix constructed from 
this sample. A procedure is developed to correct 
for this effect. 

Taking explicitely into account the number k 
of fake photons in a r decay, the efficiency matrix 
is taken as £jik , i. e. the efficiency for a produced 
event in class j with k fake photons reconstructed 
in class i. In this way the effect of the different 
fake photon multiplicities in data and simulation 
can be explicitly taken into account, assuming 
that fake photons are randomly produced. Then, 
it is sufficient to determine one factor f 3 in each 
produced class j in order to quantify the deficit in 
the simulation. These factors are obtained from 
fits of the observed P 7 distributions in data and 
the number of fake photons produced in the sim¬ 
ulation. 

In the published analysis using 1991-1993 
data |j] a much less sophisticated approach was 
taken. The deficit of fake photons in the sim¬ 
ulation was determined in a global way to be 
(16±8)%, common to all channels. Since it would 
have been delicate to generate extra fake photons, 
the procedure chosen was to actually do the op¬ 
posite, i.e. randomly kill identified (in the sense 
of the matching procedure discussed above) fake 
photons in the simulation, determine the new ef¬ 
ficiency matrix and take as corrections the sign- 
reversed deviations. 

It is clear that the current way of dealing with 
the fake photon problem is both more precise and 
more reliable. The fact that different channels 
are treated separately provides a handle on the 
different origins of the fake photons, since, for 
example, fake photons in the h class only orig¬ 
inate from hadronic interactions, whereas they 
come from both hadronic interactions and pho¬ 
ton shower fluctuations in the hn° class. 

3. Results 


2.8. Adjusting the number of fake photons 
in the simulation 

As could be expected, the number of fake pho¬ 
tons in the simulation does not agree well with 
the rate in data. Indeed, it is observed that the 
rate of low-energy simulated fake photons is insuf¬ 
ficient. This effect will affect the classification of 
the reconstructed final states in the Monte Carlo 


3.1. Determination of the branching ratios 

The branching ratios are determined using 


n° bs -n* k9 ~ 


= 5>a'T‘ 


Bj = 


lvrProd 

1 J 

Ej NY od 


(1) 

( 2 ) 
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where n° bs is the observed events number in re¬ 
constructed class i, n bkg the non-r background in 
reconstructed class i. £ji the efficiency of a pro¬ 
duced class j event reconstructed as class i. and 
jyprod p roc juced events number of class j. The 
class numbers i, j run from 1 to 14, the last one 
corresponding to the rejected r candidates. 

The efficiency matrix c ?7 ; is determined from the 
Monte Carlo, but corrected using data for many 
effects such as particle identification efficiency, rr 
selection efficiency and fake photon simulation. 

The analysis assumes a standard r decay de¬ 
scription. One could imagine unknown decay 
modes not included in the simulation, but since 
large detection efficiencies are achieved in the rr 
selection which is therefore robust, these decays 
would be difficult to pass unnoticed. An indepen¬ 
dent measurement of the branching ratio for un¬ 
detected decay modes, using a direct search with 
a one-sided r tag, was done in ALEPH 0 , limit¬ 
ing this branching ratio to less than 0.11% at 95% 
CL. This result justifies the assumption that the 
sum of the branching ratios for visible r decays is 
equal to one. 

The branching ratios are obtained and listed in 
Table g. 

3.2. Determination of systematic uncer¬ 
tainties 

Wherever possible the efficiencies relevant to 
the analysis have been determined using ALEPH 
data, either directly on the rr sample itself or on 
specifically selected control samples, as for exam¬ 
ple in the case of particle identification. The re¬ 
sulting efficiencies are thus measured with known 
statistical errors. 

In some cases the procedure is less straight¬ 
forward and involves a model for the systematic 
effect to be evaluated. An important example is 
given by the systematics in the simulation of fake 
photons in ECAL. In such cases the evaluation of 
the systematic error not only takes into account 
the statistical aspect, but also some estimate of 
the systematics involved in the assumed model. 
The latter is obtained from studies where the rel¬ 
evant parameters are varied in a range consistent 
with the comparison between data and Monte 
Carlo distributions. 


Quite often a specific cut on a given variable 
can be directly studied. The comparison between 
the corresponding distributions, respectively in 
data and Monte Carlo, provides an estimate of 
a possible discrepancy whose effect would be to 
change the assumed efficiency of the cut. If a 
significant deviation is observed, the correction 
is applied to the simulation to obtain the nom¬ 
inal branching ratio results, while the error on 
the deviation provides the input to the evalua¬ 
tion of the systematic uncertainty. The analysis 
is therefore repeated with a full re-classification 
of all the measured r decay candidates, changing 
the incriminated variable by one standard devi¬ 
ation. Since the new samples slighly differ from 
the nominal ones because of feedthrough between 
the different channels, the modified results are af¬ 
fected both by the systematic change in the vari¬ 
able value and the statistical fluctuation from the 
event sample which is uncommon to both sam¬ 
ples. In this case the final systematic uncertainty 
is obtained by adding linearly the modulus of the 
systematic deviation observed and the statisti¬ 
cal error from the monitored uncommon sample. 
This procedure is followed for all the systematic 
studies. 

Finally, the systematic deviations for each 
study are obtained with their sign in each mea¬ 
sured decay channel, thus providing the full in¬ 
formation on the correlations between the results 
and allowing the corresponding covariance matrix 
to be constructed. 

The most important systematic uncertainties 
originate from photon identification and n° re¬ 
construction, and secondary interactions. 

3.3. From reconstructed classes to exclu¬ 
sive modes 

So far branching fractions have been deter¬ 
mined in 13 classes corresponding to major r 
decay modes. However, these classes still con¬ 
tain the contributions from final states involv¬ 
ing kaons. The latter are coming from Cabibbo- 
suppressed r decays or modes with a KK pair, 
both characterized by small branching ratios com¬ 
pared to the nonstrange modes without kaons. 
Complete analyses of r decays involving neutral 
or charged kaons have been performed by ALEPH 


Table 2 

Branching ratios (%) for the reconstructed topologies (quasi-exclusive modes) from 1991-1993 and 1994- 


1995 data sets; the first error is statistical and the second is systematic. 


class 

91-93 

94-95 

e 

17.859 ± 0.112 ± 0.058 

17.799 ± 0.093 ± 0.045 

v 

17.356 ± 0.107 ± 0.055 

17.273 ± 0.087 ± 0.039 

h 

12.238 ± 0.105 ± 0.104 

12.058 ± 0.088 ± 0.083 

h,TT° 

26.132 ± 0.150 ± 0.104 

26.325 ± 0.123 ± 0.090 

h2iT 0 

9.680 ± 0.139 ± 0.124 

9.663 ± 0.107 ± 0.105 

h3ir° 

1.128 ± 0.110 ± 0.086 

1.229 ± 0.089 ± 0.068 

Mi r° 

0.227 ± 0.056 ± 0.047 

0.163 ± 0.050 ± 0.040 

3 h 

9.931 ± 0.097 ± 0.072 

9.769 ± 0.080 ± 0.059 

3h 7T° 

4.777 ± 0.093 ± 0.074 

4.965 ± 0.077 ± 0.066 

3h2n° 

0.517 ± 0.063 ± 0.050 

0.551 ± 0.050 ± 0.038 

3h3n 0 

0.016 ± 0.029 ± 0.020 

-0.021 ± 0.023 ± 0.019 

5 h 

0.098 ± 0.014 ± 0.006 

0.098 ± 0.011 ± 0.004 

5/l7T° 

0.022 ± 0.010 ± 0.009 

0.028 ± 0.008 ± 0.007 

14 

0.017 ± 0.043 ± 0.042 

0.099 ± 0.035 ± 0.037 


on the full LEP I data |||i]||. They are sum¬ 
marized and measurements with Kg and Ai£ are 
combined in Ref. |)|]. 

The r decays involving 77 or to mesons also re¬ 
quire special attention in this analysis because of 
their electromagnetic decay modes. Indeed the 
final state classification relies in part on the ir° 
multiplicity, thereby assuming that all photons — 
but those specifically identified as bremsstrahlung 
or radiative— originate from 7 r° decays. There¬ 
fore the non-7r° photons from rj and lo decays 
are treated as n° candidates in the analysis and 
the systematic bias introduced by this effect must 
be evaluated. The corrections are based on spe¬ 
cific measurements by ALEPH of r decay modes 
containing those mesons [|h|. Thus the final re¬ 
sults correspond to exclusive branching ratios ob¬ 
tained from the values measured in the topolog¬ 
ical classification, corrected by the removed con¬ 
tributions from AT, 77 and cu modes measured sep¬ 
arately, taking into account through the Monte 
Carlo their specific selection and reconstruction 
efficiencies to enter the classification. This book¬ 
keeping takes into account all the major decay 
modes of the considered mesons including 
the isospin-violating w —> 7 r + 7 r _ decay mode. 
The main decay modes into consideration are 


7tu, and 7 r 7 r °?7 with branching fractions of 

(2.26 ±0.18) 10~ 2 , (4.3 ±0.5) 10" 3 , and (1.80 ± 
0.45) 10 -3 (l 8 ), respectively. The first two num¬ 
bers are derived from the branching ratios for the 
37 T 7 T 0 and 37r27r° modes obtained in this analysis 
and the measured u> fractions of 0.431±0.033 from 
ALEPH @ and the average value, 0.78 ± 0.06, 
from ALEPH Q and CLEO Q, respectively. 

Some much smaller contributions with 77 and u> 
have been identified and measured by CLEO JhJ 
with the decay modes r —> r/n~Tr + n~ ((2.4 ± 
0.5) 10" 4 ), t r]TT~2ir° ((1.5 ± 0.5) lO’ 4 ), r -4 
W7r - 7T + 7T - ((1-2 ± 0.2) 10 -4 ), and r —> W7r _ 27r° 
((1.5 ± 0.5) 10~ 4 ). Even though the corrections 
from these channels are very small they have been 
included for the sake of completeness. Finally, 
another very small correction has been applied 
to take into account the ai radiative decay into 
TT'y with a branching fraction of (2.1 ± 0.8) 10 -3 
obtained from Ref. ( 22 ). 

3.4. Overall consistency test 

Rejected r hemispheres because of charged par¬ 
ticle identification cuts (2 GeV minimum mo¬ 
mentum and ECAL crack veto for some 1-prong 
modes, strict definition of higher multiplicity 
channels) are placed in class 14. As already em- 
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phasized, this sample does not correspond to a 
nominal r decay mode and should be explained by 
all other measured fractions in the other classes 
and the efficiency matrix. Thus the determina¬ 
tion of a hypothetical signal in this class is a 
measure of the level of consistency achieved in 
the analysis. 

For this determination the efficiency of the 
possible signal in class 14 is taken to be 100%. 
The results, already shown in Table 0 sepa¬ 
rately for the 1991-1993 and 1994-1995 data 
sets, are consistent and are combined to give 
B u = (0.066 ± (0.027) 

stat ± (0.021) 8 „ 8tiC ± 

(0.025 ) sys t,unc)%, where the last two errors re¬ 
fer to the common and uncommon uncertainties 
from the two data sets. With a combined er¬ 
ror of 0.042% this value is consistent with zero 
and provides a good check of the overall proce¬ 
dure at the nontrivial 0 . 1 % level for branching 
ratios. It is interesting to note that this value 
coincides, approximately and accidentally, with 
the limit achieved of 0.11% at 95% CL in a di¬ 
rect search for ’invisible’ decays not selected in 
the 13-channel classification. 

In the following it is assumed that all r decay 
modes have been properly considered at the 0 . 1 % 
precision level and no physics contribution be¬ 
yond standard r decays is further allowed. Thus 
the quantity B 14 is now constrained to be zero. 

It can be further noticed that this analysis pro¬ 
vides a branching ratio in the 37 r 37 r° class which 
is consistent with zero for both 1991-1993 and 
1994-1995 data sets (see Table 0 ). The result is 
therefore given as an upper limit at 95% CL 

B 3v3 *0 < 4.9 10 - 4 (3) 

consistent with the measurement made by 
CLEO H yielding B 3n3n 0 = (2.2 ± 0.5) 10" 4 . 
The final state is dominated by rj and u> reso¬ 
nances Hj| and using other channels allows one 
to obtain a lower limit for this branching ratio, 
(2.6 ± 0.4) 10" 4 . In the following a value of 
(3 ± 1) 10 -4 is used as input for this channel and 
the global analysis is performed in terms of the 
remaining 12 defined channels which are refitted. 

Figs. 0 i I an d 0 present a few comparisons 
between data and simulation, providing global 
checks on the data quality. 



p(GeV/c) 


Figure 4. Momentum spectrum in the electron 
channel in data (error bars) and Monte Carlo 
(histogram) for the 1994-1995 data set. The 
shaded histogram is the contribution of non-r 
background. 


3.5. Comparison of 1991-1993 and 1994- 
1995 results 

Since the same procedure is applied for the 
analyses of 1991-1993 and 1994-1995 data the re¬ 
sults must be consistent within the statistical er¬ 
rors of data and Monte Carlo. Good agreement 
is observed with a x 2 of 8.6 for 12 DF. In conclu¬ 
sion, the two independent data and Monte Carlo 
samples give consistent results. 
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p(GeV/c) 


Figure 5. Momentum spectrum in the muon 
channel in data (error bars) and Monte Carlo 
(histogram) for the 1994-1995 data set. The 
shaded histogram is the contribution of non-r 
background. 


3.6. Final combined results 

Finally the two sets of results are combined. 
Whether one uses only statistical or total weights 
—in the latter taking into account correlated er¬ 
rors from dynamics and secondary interactions— 
gives almost identical results. Using the total 
weights one obtains the final results shown in Ta¬ 
ble j|. 

The branching ratios obtained for the different 
channels are correlated with each other. On one 
hand the statistical fluctuations in the data and 
the Monte Carlo sample are driven by the rnulti- 




Ey (GeV) 


Figure 6. Comparison of photon energy spectra 
for all r events in data (error bars) and Monte 
Carlo (histogram) after all the corrections for 
1994-1995 data set. The dotted line shows the 
Monte Carlo distribution before correction and 
the shaded histogram is the contribution of non- 
r background. 


nomial distribution of the corresponding events, 
producing well-understood correlations . On the 
other hand the systematic effects also induce im¬ 
portant and nontrivial correlations between the 
different channels. All the systematic studies 
were done keeping track of the correlated vari¬ 
ations in the final branching ratio results, thus 
allowing a proper propagation of errors. 

The present results are consistent with those 
of the previously published ALEPH analysis [ll 
j[|. The leptonic branching ratios also agree 
within errors with the results of an independent 
ALEPH analysis which does not rely on a global 
method |24jj . 
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Invariant mass (GeV/c 2 ) 


Figure 7. Comparison of mass spectra in data 
(error bars) and Monte Carlo (histogram) after 
fake photon correction for the ft, 7 r°sample in 1994- 
1995 data sample. The dark shaded histogram 
is the contribution of non-r background and the 
light shaded is from r feedthrough. 


4. Discussion of the results 

4.1. Comparison with other experiments 

Figs, g, g, 

analysis are in good agreement with those from 
other most precise experiments. In all these cases, 
ALEPH achieves the best precision. 

A meaningful comparison can be performed be¬ 
tween the exclusive fractions and the topological 
branching ratios. Even though the latter have 
essentially no physics interest, their determina¬ 
tion can constitute a valuable cross check as they 
depend only on selection efficiency, tracking, han¬ 
dling of secondary interactions and electron iden¬ 
tification for photon conversions, and not on pho¬ 


Table 3 

Combined results for the exclusive branching ra¬ 
tios (B) for modes without kaons. The contribu¬ 
tions from channels with 77 and u> are given sep¬ 
arately, the latter only for the electromagnetic ui 
decays. All results are from this analysis, unless 
explicitly stated. The values labelled * and **,*** 
are taken from ALEPH [(Tb) and CLEO pl| , ^o| , 
respectively. 


mode 

B ^stat ^ ^syst [%] 

e 

17.837 

± 

0.072 ± 0.036 


17.319 

± 

0.070 ± 0.032 

TV~ 

10.828 

± 

0.070 ± 0.078 

7r _ 7r° 

25.471 

± 

0.097 ± 0.085 

3 

1 

to 

0 

9.239 

± 

0.086 ± 0.090 

0 

CO 

1 

0.977 

± 

0.069 ± 0.058 

0 

1 

0.112 

± 

0.037 ± 0.035 

+ 

1 

1 

9.041 

± 

0.060 ± 0.076 

0 

+ 

1 

1 

4.590 

± 

0.057 ± 0.064 

7T 7T 7T + 27r° 

0.392 

± 

0.030 ± 0.035 

7 r _ 7 r _ 7 r + 3 - 7 r° (estim.) 

0.013 

± 

0.000 ± 0.010 

CO 

1 

to 

+ 

0.072 

± 

0.009 ± 0.012 

CO 

1 

to 

+ 

0 

0.014 

± 

0.007 ± 0.006 

— 0 * 

7T 7T Ty 

0.180 

± 

0.040 ± 0.020 

— 0 0 * * 

n 2n rj 

0.015 

± 

0.004 ± 0.003 

— — 4 - ** 
tv tv tv ij 

0.024 

± 

0.003 ± 0.004 

a 3 (—> 7 r - 7 ) (estim.) 

0.040 

± 

0.000 ± 0.020 

> 7T°7,7T + 7r _ )* 

0.253 

± 

0.005 ± 0.017 

TV TV LU{—> TV 7 , 7T TV ) ’ 

0.048 

± 

0.006 ± 0.007 

— 0 0 / 0 4- — \ ** 

tv 2tv —> tv 7 , tv tv ) 

0.002 

± 

0.001 ± 0.001 

TV TV TV COy—> TV 7 , TV TV ) 

0.001 

± 

0.001 ± 0.001 


ton identification. The results from this analysis 
can be compared in this way with a dedicated 
analysis recently performed by DELPHI p5| . 
Summing up appropriately (both analyses assume 
a negligible contribution from hadronic multiplic¬ 
ities higher than 5, in agreement with the 90% CL 
limit by CLEO J2(j] one gets Bj < 2.4 10 -6 ) 

B 3 = (14.652 ± 0.067 ± 0.086) % (4) 

B 5 = (0.093 ± 0.009 ± 0.012) % (5) 

in good agreement with the DELPHI values, 
B 3 = (14.569 ± 0.093 ± 0.048) % and B 5 = 
(0.115 ± 0.013 ± 0.006) %. The rather small sys¬ 
tematic uncertainty in the DELPHI results re¬ 
flects the fact that a sharper study of hadronic 


10 . 11 show that the results of this 
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-ALEPH 89-90 

18 . 09 ± 0 . 45 ± 0.45 

ARGUS 

17 . 5 ± 0 . 3 ± 0.5 

ALEPH 91-93 

17 . 79 ± 0 . 12 ± 0.06 

CLEO 97 

17 . 76 ± 0 . 06 ± 0.17 

DELPHI 91-95 

17 . 877 ± 0 . 10910.11 

OPAL 91-95 

17 . 8110 . 0910.06 

L3 91-95 

17 . 80610 . 10410.076 

ALEPH 91-95 Prel. 

..17.837+Q.Q72+Q.Q36. 

average 

17 . 823 + 0.051 


18 19 

B(evv) (%) 
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ALEPH 89-90 

17 . 3510 . 4110.37 

ARGUS 

17 . 410 . 310.5 

OPAL 90-92 

17 . 3610.2710 

ALEPH 91-93 

17 . 3110 . 1110.05 

CLEO 

17 . 3710 . 0810.18 

DELPHI 91-95 

17 . 32510 . 09510.077 

L3 91-95 

17 . 34210 . 1110.067 

ALEPH 91-95 Prel. 

. 1 . 7 . 31 . 9 + 0 . 07 + 0.032 . 

average 

17 . 331 + 0.054 


17 18 19 

B(|ivv) (%) 


Figure 8 . Comparison of ALEPH measurement 
with published precise results from other experi¬ 
ments for t —► evD. 


Figure 9. Comparison of ALEPH measurement 
with published precise results from other experi¬ 
ments for t —> fjivv. 


interactions can be performed when only charged 
particles are considered in the analysis. In ad¬ 
dition the modes with K® —+ 7 r + 7 r _ decays are 
subtracted statistically here, rather than trying 
to identify them on an event-by-event basis. 


4.2. Universality in the leptonic charged 
current 

4.2.1. fi-e universality from the leptonic 
branching ratios 

In the standard V-A theory with leptonic cou¬ 
pling gi at the WWj vertex, the r leptonic par¬ 
tial width can be computed, including radiative 
corrections ]27j and safely neglecting neutrino 
masses: 


r(r -> v T h>i( 7 )) 


G T Gim 5 T f mf 
1927T 3 \m% 


SWZ, 


( 6 ) 


where 

Gi = 


9i 


4V2 


&w 

s; 

fix) 


1 + 


3 ml 


1 + 


5 Myy 
x(m T ) 
2tt 



1 — 8x + 8 x 3 — x 4 — 12 a; 2 ln.T 


( 7 ) 


Numerically, the W propagator correction and 
the radiative corrections are small: = 1 + 

2.9 • 1CT 4 and 6* = 1 + 43.2 • 1CT 4 . 

Taking the ratio of the two leptonic branching 
fractions, a direct test of g-e universality is ob¬ 
tained. The measured ratio 

D 

-Af = 0.9709 ± 0.0060 ± 0.0029 (8) 

Be 

agrees with the expectation which is equal to 
0.97257 when universality holds. Alternatively 
the measurements yield the ratio of couplings 

^ = 0.9991 ±0.0033 (9) 

9e 

which is consistent with one. 
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B(hv) (%) B(h^°v) (%) 


Figure 10. Comparison of ALEPH measurement 
with published precise results from other experi¬ 
ments for r —> hu (sum of iris and Ku). 


Figure 11. Comparison of ALEPH measurement 
with published precise results from other experi¬ 
ments for t — > hn°v (sum of itk 0 v and Kn°u). 


This result is in agreement with the best test 
of /u-e universality of the W couplings obtained 
in the comparison of the rates for n —> and 

7 r —> ev e decays, where the results of the two most 
accurate experiments j2^,^9| can be averaged to 
yield ^ = 1.0012 ±0.0016. The results have 
comparable precision, but it should be pointed 
out that they are in fact complementary. The r 
result given here probes the coupling to a trans¬ 
verse W (hclicity ±1) while the 7 r decays measure 
the coupling to a longitudinal W (helicity 0). It 
is indeed conceivable that either approach could 
be sensitive to different nonstandard corrections 
to universality. 

Since B e and B^ are consistent with /i-e univer¬ 
sality their values can be combined, taking com¬ 
mon errors into account, into a consistent leptonic 

branching ratio for a massless lepton (essentially 

2 

the case for the electron, noting that /(^f) = 1 


for all practical purposes) 

B ( mi =o) = 17 _ 822 ± o.044 ± 0.022(%) , (10) 

where the first error is statistical and the second 
from systematic effects. 


4.2.2. Tests of t-/j, and r-e universality 

Comparing the rates for T(r —> is T eV e ( 7 )), 
P(r —*• 7 )) and T(n —> v ll eV e ('y)) provides 

direct tests of the universality of r-/i-e couplings. 
Taking the relevant ratios with calculated radia¬ 
tive corrections, one obtains 







<1 

<1 

(11) 

■ Aw 

(12) 
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where /(- h = 0.9998, A w = = 1-2.9-10" 4 , 

A 7 = = 1 + 8.5 • 10 -5 , and t ; is the lepton l 

lifetime. 

From the present measurements of f? e , B the 
r mass (Ig| ], m r = (1777.031 q; 26) MeV (domi¬ 
nated by the BES result @), the r lifetime Jl9| , 
r r = (290.6± 1.1) fs and the other quantities from 
Ref. |ll|, universality can be tested: 

— = 1.0009 ± 0.0023 ±0.0019 ±0.0004 (13) 

— = 1.0001 ± 0.0022 ±0.0019 ±0.0004 ,(14) 

9e 

where the errors are from the corresponding lep- 
tonic branching ratio and the r lifetime and mass, 
respectively. 


4.2.3. r-/i universality from the pionic 
branching ratio 

The measurement of B w also permits an inde¬ 
pendent test of t-/j, universality through the rela¬ 
tion 


B* r n 2 m^ml 


B 


-JT— t t TTI t 

1 - ml/ml 
1 — ml / m1 


/7T 5 


(15) 


where the radiative correction |3l| amounts to 
S T /n = 1.0016 ± 0.0014 Using the world-averaged 
values for the r and 7r (t t and ry) lifetimes, and 
the branching ratio for the decay n —> /j,v ]l9| , 
the present result for B n , one obtains 


— = 0.9962 ± 0.0048 ± 0.0019 ± 0.0002 , (16) 

9» 

comparing the measured value (B , T = 10.823 ± 
0.104)% to the expected one assuming universal¬ 
ity (10.910±0.042)%. The quoted errors in Eq. |l6| 
are from the pion mode branching ratio and the 
r lifetime and mass, respectively. 

The two determinations of — obtained from B e 
and Bn are consistent with each other and can be 
combined to yield 

— = 1.0000 ± 0.0021 ± 0.0019 ± 0.0004 , (17) 

9n 


ALEPH Preliminary 


' » ’ 1 » ’ 


1 1 1 1 1 1 1 1 

- 1 B e ALEPH 



— B u ALEPH 



B ALEPH 

7t 



-h x T WA 

_1_1_1_1_1_1_ 

_l_ 

average 

_._._1_I_1_,_1_■_ 


17.4 17.6 17.8 18 18.2 


B e (%) 


Figure 12. The measured value for B e compared 
to predictions from other measurements assuming 
leptonic universality. The vertical band gives the 
average of all determinations. 


where the errors are from the electron and pion 
branching ratio and the r lifetime and mass, re¬ 
spectively. Universality of the r and fj, charged- 
current couplings holds at the 0.29% level with 
about equal contributions from the present deter¬ 
mination of B e and and the world-averaged 
value for the r lifetime. 

The consistency of the present branching ratio 
measurements with leptonic universality is dis¬ 
played in Fig. |l2| where the result for B e is com¬ 
pared to computed values of B e using as input B M 
(assuming e — /i universality), ry and (fi — t 
universality), and B w and r 7r (/x — r universality). 
All values are consistent and yield the average 

^universality = ( 1? 81Q ± Q.039) % . (18) 

4.2.4. decays to 37r and 7 t27t 0 

With the level of precision reached it is interest¬ 
ing to compare the rates in the 37 t and 7r27r° chan- 
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nels which are completely dominated by the a± 
resonance. The dominant pir intermediate state 
leads to equal rates, but a small isospin-breaking 
effect is expected from different charged and neu¬ 
tral 7 T masses, slightly favouring the tt2tt° chan¬ 
nel. 

A recent CLEO partial-wave analysis of the 
final state |l 2 j has shown that the sit¬ 
uation is in fact much more complicated with 
many intermediate states, in particular involv¬ 
ing isoscalars, amounting to about 20 % of the 
total rate and producing strong interference ef¬ 
fects. A good description of the a\ decays was 
achieved in the CLEO study, which can be ap¬ 
plied to the 37 r final state, predicting j3^] a ratio 
of the rates 3tt/tt2i r° equal to 0.985. This value, 
which includes known isospin-breaking from the 
pion masses, turns out to be in good agreement 
with the measured value from this analysis which 
shows the expected trend 

Bo 

— = 0.979 ± 0.018 . (19) 

J-'Tv2n° 

4.2.5. The tttt 0 branching ratio in the con¬ 
text of a t }, ad 

The 7 T 7 T 0 final state is dominated by the p res¬ 
onance as demonstrated in Fig. [f]. Its mass dis¬ 
tribution —or better the corresponding spectral 
function— is a basic ingredient of vacuum polar¬ 
ization calculations, such as that used for com¬ 
puting the hadronic contribution to the anoma¬ 
lous magnetic moment of the muon a dad . In this 
case the p contribution is dominant (71%) and 
therefore controls the final precision of the result. 
The current evaluation j33| used by the BNL ex¬ 
periment [[ii] is based on r and e + e _ data, with 
more precision on the r side. 

The normalization of the spectral function is 
provided by the branching fraction B^ n o taken 
from the world average and completely domi¬ 
nated by the published ALEPH result The 
new result given here is larger by 0 . 68 %, thus one 
can expect a slightly larger contribution to a ^, ad , 
but it remains within the quoted uncertainty ||i.'i| 
of 6.2 • 10 -10 , corresponding to a relative error of 
1 . 2 % for the p contribution. 

A new evaluation |35|,[i6| is available, using 
the spectral functions from the present prelimi- 



B -B (°/' 

“measured “predicted ' • 


Figure 13. The differences between the measured 
values for and B ^ 27r o and their respective pre¬ 
dictions from leptonic universality and isospin 
symmetry with 7r-mass breaking, compared to the 
precision on B nv o. 


nary ALEPH analysis, the published CLEO re¬ 
sults Q and new results from e + e“ annihila¬ 
tion from CMD-2 |Q. Since a disagreement is 
observed between the r and e + e“ spectral func¬ 
tions, it is important to check all ingredients, in 
particular the determination of the branching ra¬ 
tio B nv o. As most of the systematic uncertainty 
in B nir o comes from 7 / 71 "° reconstruction, it is in¬ 
teresting to cross check the results in the ’adja¬ 
cent’ hadronic modes, i.e. the n and 7 t27t 0 chan¬ 
nels. This is possible if universality in the weak 
charged current is assumed, leading to an abso¬ 
lute prediction o f B^ using as input the r life¬ 
time (see Section 1.2.2| ), and by computing B n2 ti - 0 
from the measurement of B^ which is essentially 
uncorrelated with o (see Section 4.2.4 ). The 
comparison, shown in Fig. [l3|, does not point to 
any systematic bias in the determination of B^^o 
within the quoted uncertainty. 
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4.2.6. Separation of vector and axial- 
vector contributions 

From the complete analysis of the r branching 
ratios presented in this paper, it is possible to 
determine the nonstrange vector (V) and axial- 
vector (A) contributions to the total r hadronic 
width, conveniently expressed in terms of their 
ratios to the leptonic width, called R T y and R t ,a, 
respectively. The determination of the strange 
counterpart R T y is already published |J. 

The ratio R T for the total hadronic width is 
calculated from the difference of the ratio of the 
total hadronic width and electronic branching ra¬ 
tio, 

= 1 -Be-Bp = 1 _ 1 97257 

B e B e 

= 3.642 ±0.012. (20) 

taking for B{r~ —> e~ i' e r , T ) the value obtained in 
Section |4.2| assuming universality in the leptonic 
weak current. Using the ALEPH measurement of 
the strange width ratio ||, very slightly modified 
to take into account the channel K*~ij measured 
by CLEO H 

R t S = 0.1625 ±0.0066 , (21) 

the following result is obtained for the nonstrange 
component 

R T y+A = 3.480 ± 0.014 . (22) 

Separation of V and A components in hadronic 
final states with only pions is straightforward. 
Isospin invariance relates the spin-parity of such 
systems to their number of pions: G-parity =1 
(even number) corresponds to vector states, while 
G=-l (odd number) tags axial-vector states. This 
property places a strong requirement on the effi¬ 
ciency of 7r° reconstruction, a constraint that was 
strongly integrated in this analysis. 

Modes with a KK pair are not in general eigen¬ 
states of G-parity and contribute to both V and A 
channels. The respective components have been 
determined in the ALEPH analysis |)jj. While the 
decay to I\~K° is pure vector, the KKtt mode 
has been shown to be almost completely axial vec¬ 
tor, with a fraction 0.941°'^. This information 
was not available at the time of the previous anal¬ 


ysis of the nonstrange modes [|| where a conser¬ 
vative value of 0.5 ± 0.5 was used. For the de¬ 
cays into KKirir no information is available in 
this respect and the same conservative fraction is 
assumed. 

The total nonstrange vector and axial-vector 
contributions obtained in this analysis are: 

Rr V = 1.778 ± 0.010 ± 0.002 , (23) 

R r A = 1.701 ± 0.011 ± 0.002 , (24) 

where the second errors reflect the uncertainties 
in the V/A separation in the channels with KK 
pairs. Taking care of the correlations between the 
respective uncertainties, one obtains the differ¬ 
ence between the vector and axial-vector compo¬ 
nents, which is physically related to the amount of 
nonperturbative QCD contributions in the non¬ 
strange hadronic r decay width: 

R T y_ A = 0.077 ±0.018 ±0.005 , (25) 

where again the second error has the same mean¬ 
ing as in Eqs. ( f±i[ ) and (|24|). The ratio 

Br.y-A = o Q22 ± 0.005 (26) 

Rt,V+A 

is a measure of the relative importance of non¬ 
perturbative QCD contributions. 

5. Conclusions 

A final analysis of r decay branching fractions 
using all LEP-I data with the ALEPH detector 
is presented. As in the publication based on 
the 1991-1993 data it uses a global analysis of 
all modes, classified according to charged par¬ 
ticle identification, and charged particle and tt° 
multiplicity up to 4 7r°s in the final state. Ma¬ 
jor improvements are introduced with respect to 
the published analysis and a better understand¬ 
ing is achieved, in particular in the separation 
between genuine and fake photons. In this pro¬ 
cess shortcomings and small biases of the previ¬ 
ous method were discovered are corrected, lead¬ 
ing to more robust results. As modes with kaons 
(A' ± , Kg, and Kg) have already been studied and 
published with the full statistics, the nonstrange 
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Table 4 

A summary list of ALEPH branching ratios (%). The labels V, A and S refer to the nonstrange vector 
and axialvector, and strange components, respectively. 


mode 

B ±o- to t [%] 

(PRELIMINARY!) 


e 

17.837 ± 0.080 




17.319 ± 0.077 



7T~ 

10.828 ± 0.105 

A 


TT~TT° 

25.471 ± 0.129 

V 


o 

CM 

1 

9.239 ± 0.124 

A 


7 r 37 t° 

0.977 ± 0.090 

V 


O 

1 

0.112 ± 0.051 

A 


+ 

9.041 ± 0.097 

A 


+ 

o 

4.590 ± 0.086 

V 


7r 7T 7T + 27r° 

0.392 ± 0.046 

A 


+ 

co 

o 

0.013 ± 0.010 

V 

estim 

37r _ 27r + 

0.072 ± 0.015 

A 


37r _ 27r + 7r° 

0.014 ± 0.009 

V 


TT~TT°TI 

0.180 ± 0.045 

V 


(3tt)-77 

0.039 ± 0.007 

A 

CLEO 


0.040 ± 0.020 

A 

estim 

7T _ Uj (—> 7r 0 7,7T + 7r _ ) 

0.253 ± 0.018 

V 


7r _ 7r°a;(—» 7T°7,7r + 7r~) 

0.048 ± 0.009 

A 

+ CLEO 

(37t) — cu(—* 7r°7,7r + 7r _ ) 

0.003 ± 0.003 

V 

CLEO 


0.163 ± 0.027 

V 


K~tt°K° 

0.145 ± 0.027 

(94±<>)%A 


n~K°K° 

0.153 ± 0.035 

(94+®)%A 


K-K+tt- 

0.163 ± 0.027 

(94±®)%A 


(KKtttt)- 

0.05 ± 0.02 

(50 ± 50)% A 


K~ 

0.696 ± 0.029 

S 


K~tt° 

0.444 ± 0.035 

S 


K°n~ 

0.917 ± 0.052 

S 


K~ 27t° 

0.056 ± 0.025 

S 


I\ _ 7r + 7T _ 

0.214 ± 0.047 

S 


A' 0 7T _ 7r 0 

0.327 ± 0.051 

S 


(K3tt)- 

0.076 ± 0.044 

S 


K~r) 

0.029 ± 0.014 

S 
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modes without kaons are emphasized. Taken to¬ 
gether these results provide a complete descrip¬ 
tion of r decay modes up to 6 hadrons in the 
final state. 

The measured branching ratio values are in¬ 
ternally consistent and agree with known con¬ 
straints from other measurements in the frame¬ 
work of the Standard Model (or even looser 
assumptions). The precision reached and the 
completeness of the results are for the moment 
unique. More specifically, the results on the lep- 
tonic and pionic fractions lead to powerful tests 
of universality in the charged leptonic weak cur¬ 
rent, showing that the e — /i — r couplings are 
equal within 2-3 per mille. The branching ra¬ 
tio of r —> v t txtC which is of particular inter¬ 
est to the accurate determination of vacuum po¬ 
larization effects is determined with a precision 
of 0.5% to be (25.47 ± 0.13) %. Also the ratio 
of eti branching fractions into n2n° and 37 t final 
states is measured to be 0.979 ± 0.018, in agree¬ 
ment with expectation from partial wave analyses 
of these decays. Separating nonstrange hadronic 
channels into vector (V) and axial-vector (A) 
components and normalizing to the electronic 
width yields the ratios R T y = 1.778 ± 0.010, 
R t ,a = 1-701 ± 0.011, R t ,v+a = 3.480 ± 0.014 
and R t ,v-a = 0.077 ± 0.019. 

The results presented here are combined with 
previously published ALEPH results on final 
states with kaons in Table [|. 
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